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ABSTRACT 


Ocean surface circulation is an essential component of the world climate system. 
In this study, the Ocean Surface Currents Analysis - Real Time (OSCAR) data, derived 
from satellite altimeter and scatterometer, is used to investigate the connection between 
the Indian Ocean Dipole and eastward equatorial jet. The raw OSCAR data set was 
refined using the Optimal Spectral Decomposition method (OSD). Data was analyzed to 
show the seasonal variability of the surface currents in the Indian Ocean with emphasis 
on the Somali current. North Arabian Sea, and Bay of Bengal. 

To investigate the link between the Eastward Equatorial jets (Wyrtki jets) and the 
Indian Ocean Dipole (lOD) Mode events, complex EOE analysis was applied to the 
currents in the Equatorial region. This analysis revealed that the spatial anomaly of the 
currents along the equator can directly relate to the occurrence of a Dipole Mode event. 
Eurther zonal currents anomalies and vector plots, for the years with Dipole Mode events, 
also depicted the anomalous behavior of the equatorial jets during a Dipole Mode event. 
This study demonstrates that satellite-based altimetry data can be used to refine the 
climatological knowledge of a certain region in the world. 
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I. 


INTRODUCTION 


A. GENERAL INFORMATION 

Much work has been done in the past to understand the surface circulation in the 
Indian Ocean. One of the most interesting features about the Indian Ocean climatology is 
its seasonal reversal or, in other words, the monsoon cyele. Its surface circulation is 
highly analogous to the seasonal pattern of the winds in the region. This study uses 
satellite-based altimetry data to analyze the surfaee circulation patterns in the Indian 
Ocean region. To aceomplish this task, the observational satellite-based altimeter data 
from http://www.oscar.noaa.gov/index.html was utilized. The website contains near real¬ 
time ocean surface currents derived from satellite altimeter (JASON -1, GFO and 
ENVISAT) and scatterometer data with 1° x 1° resolution for the world oceans from 
59.5° S to 59.5° N. This data is posted online as OSCAR (Ocean Surface Current 
Analysis - Real-time) data. This is an invaluable resource as it can be used for large-seale 
climate diagnostics and prediction, fisheries management, monitoring; debris drift, larvae 
drift, oil spills, fronts and eddies. This data can also be utilized for planning search and 
rescue, naval and maritime operations such as mine warfare. The methodology of 
OSCAR combines geostrophic, Ekman and Stommel shear dynamics, and a 
complimentary term from the surface buoyancy gradient. The seeond chapter gives a 
brief synopsis of the OSCAR database. The application and uses of surface current 
analysis are discussed in the last ehapter. 

Along with all the good in the OSCAR data, there are weaknesses. One of its 
main weaknesses is the inability to reveal ocean boundary currents, such as the Gulf 
Stream and Kuroshio. They are not adequately distinguishable in the plots produced by 
the raw OSCAR data. To refine the OSCAR data, this study used the OSD (Optimal 
Spectral Decomposition) method. After reconstructing the OSCAR data through OSD, 
the resultant data produced adequately realistic surface circulation of the lateral boundary 
currents. Details about the OSD method to process the noisy OSCAR data are given in 
the third chapter. 
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One of the most dynamic research avenues, in context to Indian Ocean 
climatology, is the phenomenon of the Indian Ocean Dipole (lOD). Although occurrence 
of this climate mode can be traced to ancient history, it was not identified until 1999. This 
is a climate mode whose occurrence is dominant in the tropical regions of the Indian 
Ocean. A Dipole event brings heavy rainfall over the east African region and 
droughts/forest fires in the Indonesian region or vice versa. The fifth chapter will analyze 
certain surface circulation patterns to see how we can identify the occurrence of a dipole 
event while seeing the anomaly in surface circulation patterns along the equator. 

B. THE INDIAN OCEAN 

The Indian Ocean (Figure 1-1) is the world’s third largest ocean. It is bounded by 
Asia to the north (including the Indian subcontinent, after which it is named); by Africa 
to the west; by Indochina, Malaysia, Indonesia, and Australia to the east; and by 
Antarctica to the south. The Indian Ocean has a north-south extent of about 9600 km 
from Antarctica to the inner Bay of Bengal. It spans up to 7800 km in the east-west 
direction between Southern Africa and Western Australia. Excluding its southern part, it 
has an area of about 48x10^ km^ and, with the Southern part, the area increases up to 
74.1x10^ km^ (Tomczak & Godfrey, 1994). It separates from the Atlantic Ocean at 20° 
east meridian and from the Pacific at 147° east meridian. In the North, the Indian Ocean 
extends up to approximately 30° north latitude in the Persian Gulf.^ There are several 
small islands throughout the ocean which includes Madagascar, Comoros, Seychelles, 
Maldives, Mauritius, and Sri Lanka. 

The Indian Ocean region is the home to the world’s first known urban 
civilization: the evidence of trade through water in this region can be traced to ancient 
times. This region has remained very important in terms of its export to other parts of the 
world. In the olden times, due to its rich agriculture, its prominent product was 
agricultural goods. Although it has been the region of political interest to all the worlds’ 
strongest nations, no nation dominated this Ocean until the mid-1800s when the United 
Kingdom took control of most of its surrounding parts. Because it houses the world’s 

^ Tomczak and Godfrey, (1994). 
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most crude oil reserves, it now again has gained a lot of attention politically and 
economically. In the present time, mostly India and Australia dominate the Oeean. 
However, after the fall of the Twin Towers in U.S., the new global war against terrorism 
has resulted in US-led eoalition forees that are making a differenee in overall dominanee 
in the region. 


Indian Ocean Area 



Figure 1-1. Indian Ocean as seen on the Globe (From: www.nga.mil) 

C. POLITICAL AND ECONOMIC IMPORTANCE OF INDIAN OCEAN 

In the recent past, the Indian Ocean region has become a place of strategic 
importance to the world’s greatest powers. The reason may be that the Indian Ocean 
region has played a vital role in overall facilitation to global maritime trade. The major 
reasons whieh make Indian Ocean very important in eontext to global maritime trade are 
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crude oil reserves, population, location, historical influences, and growing economies.^ 
To see how the Indian Ocean has gained tremendous political and economic importance 
in the past few decades, there will be a brief discussion of a few major factors in the 
following subsections. 

Today the Indian Ocean houses some of the busiest waterways in the world. 
These waterways include the Suez Canal in Egypt, Bab-el-Mandeb (along Djibouti and 
Yemen border). Straits of Hormuz (along Iran and Oman border), and Straits of Malacca 
(along Indonesia and Malaysia border). All these narrow channels, or “chokepoints,” are 
critical to world crude oil trade as huge amounts of oil pass through them.3 An estimated 
40% of the world’s offshore oil production comes from the Indian Ocean.4 In political 
context, the Persian Gulf remains a crucial passageway. This is because all the world’s 
dynamic economies, which are highly dependant on foreign oil, require safe trading 
routes in and out of these passages. The dependence of the world’s biggest economies on 
foreign oil imported from the Indian Ocean rim can be judged from this fact: in present 
time, it has been estimated that European, Japanese and U.S. economies respectively 
import 70%, 76%, and 25% of their crude oil requirements from the Gulf region (Nazrey 
Khalid, 2005). China, as one of the fast- emerging and growing economies in the world, 
also has vital trading and energy interests in the Indian Ocean’s waters. The rapid 
economic progress in developing nations, such as India and South Africa, also adds up 
their interest in the safe passageways in the Indian Ocean region. ^ In short, it can be said 
that these paths not only serve as the strategic and economic lifeblood of many countries 
along them, but also play important roles as the main waterways for global trade 
especially in the energy sector. Thus, the importance of the Indian Ocean as a facilitator 
to the world’s trade cannot be undermined. 


2 Nazery Khalid (2005). 

3 World Oil Transit Chokepoints, Report by Energy Information Administration. Washington, DC. 
March 2004. 

4 World Fact Book (CIA). 

5 Nazery Khalid (2005). 
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Besides the huge oil trade in the Indian Ocean region, the vast consumer market 
makes the Indian Ocean region economically dynamic. The population of the Indian 
Ocean rim is close to two billion. This provides an enormous market for trade and 
consumerism.6 Recently there has been a significantly increasing trend: countries in this 
region moved from inward-looking economies to now opening up their economies at 
many different avenues. This development has also given rise to new opportunities for 
international trade most of which is through water. Hence, with all that has been 
discussed above about the Indian Ocean, it is evident that it will continue to play a crucial 
role in the overall facilitation to global maritime trade. 

For centuries, the Indian Ocean region has remained politically important. After 
the end of the first Gulf War in the early 1990s, its importance has increased 
exponentially in terms of providing safe passages to the maritime trade. The Indian 
Ocean rim consists of the world’s fast growing economies of India, Far East, and the Gulf 
region. All these economies, along with other big world economies, want peace and 
stability in this region to acquire safer trade. The vitality of safety of these passages can 
be judged from the consistent patrolling of U.S.-led coalition naval forces in the northern 
Indian Ocean. In short, it can be said that as the world’s bigger economies’ appetite for 
oil is growing and developing nations are coming on the road to massive 
industrialization, the Indian Ocean region will gain more and more importance 
proportionally. 

D. CLIMATOLOGY OF INDIAN OCEAN FROM PREVIOUS STUDIES 

Before we proceed to analyze the Satellite based altimeter data to investigate the 
Indian Ocean region it is necessary that we should have an idea about the general 
dynamics and climatology of this region. Amongst the general climatology it is the wind 
flow which affects the ocean surface currents the most. Therefore in the following 
subsections we will see how the conventional wind flow in the Indian Ocean region looks 
like. Most of the work presented here is taken from Regional Oceanography by Tomczak 
and Godfrey (1994). 

6 Sheth (2002). 
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One of the interesting features of the Indian Ocean wind flow pattern is its 
seasonal reversal. This is generally known as the monsoon cycle. This reversal in turn 
produces reversed climatological conditions especially in the northern Indian Ocean. The 
effects of this climatological cycle can be felt even up to southern subtropics. Figure 1-2 
summarizes the monsoon cycle of the Indian Ocean. How this seasonal reversal is 
developed, and what are its effects on the surrounding regions, is briefly discussed in this 
subsection. 
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Figure 1-2. A summary of the monsoon system in the Indian Ocean. The top part 
indicates the wind cycle; the lower part shows the major currents that 
develop in response to the wind (From: Tomczak & J. S Godfrey, 1994). 

In northern hemisphere winter (from December to March), the northeast monsoon 
determines the climatology of the northern part of the Indian Ocean. During this time, 
high pressure develops over the Indian land mass. The reason is two-fold: in winter the 
land mass cools up quickly as compared to the ocean surface and relatively hot air over 
the ocean tends to rise. This creates a zone of low pressure over the ocean. Because of 
this pressure, gradient northeasterly winds flow over the tropical region. These winds 
over the northern Indian Ocean are the Trades. However, due to their seasonal cycle they 


6 























are known as the northeast monsoon (the word ‘monsoon’ is an Arabic word meaning 
seasonally reversing winds). During this time, the Inter-tropical Convergence Zone 
(ITCZ) and the Doldrums are located south of the equator at about 5°S (Tomczak and 
Godfrey, 1994). Most of the air pressure gradient during this time is retained behind the 
Tibetan Plateau, whereas the air pressure gradient over the ocean is small. This results in 
a wind of moderate force over the ocean. This wind is also dry; therefore, the winter 
season is dry over most of southern Asia. 

During this time, in the southern hemisphere, the climatology is determined by the 
low pressure in the tropical region and by the high pressure in the sub tropical belts. As a 
result of this pressure, gradient Southeast Trades’ flow is relatively stronger and uniform 
than in the Pacific Ocean. Similar to the other oceans along the eastern coast, the wind 
flow is southerly. Alternatively, Australia’s Northwest Shelf winds become 
southwesterly. This avoids the low heat over the Australian continent. This climatological 
pattern during the southern hemisphere’s summer brings rain to northern Australia. This 
season also brings rain throughout the doldrums, with the Indonesian region getting 
maximum rainfall. 

In the northern hemisphere’s summer, it is the southwest monsoon that determines 
the climatology of the Indian Ocean. During this time of the year, the land mass heats up 
quickly due to the intense hot climate over the land compared to the ocean’s relative 
coolness. Due to this pressure dynamics, a deep low pressure develops over northern 
Arabia and Pakistan. Australia in the southern hemisphere winter now has a center of 
high pressure. The atmospheric high north of Kerguelen Islands is also shifted westward 
towards southern Africa. The pressure gradient between Arabia and Madagascar is about 
22 hPa. This never exceeds 6 hPa during the winter season (Tomczak and Godfrey, 
1994). As a result of this climatological change, the winds over the northern Indian 
Ocean completely reverse their direction and a wind jet blows in a northerly direction. 
This southwest monsoon in the northern hemisphere is the continuation of southern 
hemisphere Trades which, during this time of the year, between 10° and 20°S are 
stronger than anywhere else in the world (Tomczak and Godfrey, 1994). Over the ocean 
there is a very little change in the ITCZ but during the southwest monsoon the frequent 
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disturbances breakaway and it settles down south of the Himalayas. The southwest 
monsoon passes over the low in the Indian sub-continent to bring the well-known 
monsoon rains and floods to the region, which is very critical to the agriculture of this 
region. 
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II. OCEAN SURFACE CURRENTS ANALYSIS - REAL TIME 

(OSCAR) DATA BASE 


A. INTRODUCTION 

The web site dedicated to access Ocean Surface Currents Analyses - Real-time 
(OSCAR) data at http://www.oscar.noaa.gov/index.htnil is an invaluable resource that 
provides near real-time world oceans surface current fields. Subsequently, these fields 
can be used for different operational purposes. These surface current fields are 
automatically computed with the combination of gridded data derived from satellite- 
based altimeters and scatterometer (Bonjean et ah, 2004). The automatic mechanism to 
compute this data has been developed at Earth and Space Research (ESR), whereas, the 
web-based access system is developed at Pacific Marine Environmental Eaboratory 
(PMEE). Eigure 2-1 shows a snap shot of the home page of the web server. 





Near-realtime ocean surface currents 
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Eigure 2-1. The OSCAR website (Erom: www.oscar.noaa.gov) . 
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B. DATA PROCESSING 


The OSCAR near real-time surface current fields are available to the users every 
five days starting from October 1992 for the world oceans from 60 N to 60 S. Figure 2-2, 
taken from the OSCAR website, shows the flow diagram of how data is being processed 
and what organizations are participating to make it useful for the end user. The automated 
system developed at ESR collects up-to-date sea surface height (SSH), wind and sea 
surface temperature (SST) data. These are used to calculate ocean surface current fields 
that can then be accessed and downloaded from the web server. 



Operational Ocean Surface Currents from 
Satellite Surface Topography and Vector Winds 
Pilot Demonstration 
Data Flow and Processing Diagram 
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Realtime OuikScaU 

Sea Winds Data Access 

AJtmetei Eiror Corrections 




Surface Current Data Products and User Access 


Hawaii 

CoastWalch 



NESDIS 



Figure 2-2. OSCAR data processing flow diagram (From: www.oscar.noaa.gov) . 

The system gets the ocean surface topography, or the sea surface height data, 
from satellite-based altimeters JASON-1 (10-day cycle), GFO, and ENVISAT. The 
geostrophic winds’ data is derived from QUIKSCAT satellite. Initially, from 1992 to 
2006, the ocean surface current fields were derived from satellite-based altimeters 
TOPEX/POSEIDON and ERSl-2 and radiometer SSM/I. Downloadable data files from 
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the website ean be used for many subsequent seientifie and operational purposes. This 
may inelude seareh and reseue planning, mine warfare elimate diagnosties, and 
monitoring drifts of different suspended partieles. In many studies, this data has been 
used to monitor large-seale eddies. 

C. THE INTERACTIVE WEB-BASED ACCESS 

The interaetive web-based aeeess to OSCAR data uses JavaSeript-based 
applieations. The user ean either seleet the desired region and time frame to get an online 
plot to analyze the climatological dynamics in that particular area or the data files can be 
downloaded for further analysis. The JavaScript-based user interface provides users great 
flexibility in extraction or display of the data, as their needs demand. Users can select 
latitude-longitude maps, latitude-time sections, time-longitude sections, and time-series 
plots of desired geographical region and time period. Figures 2-3 to 2-6 show different 
user interfaces that are provided on the website. 


Data Display & Download 
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Figure 2-3. Lat-Time section plot options page (From: www.oscar.noaa.gov) . 
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Data Display & Download 


Dat* Download 



Latitude-Longitude Map for Ocean Surface Current Non-Java Version 
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Figure 2-4. Latitude-Longitude map plot options page (From: www.oscar.noaa.gov) . 


Data Display & Download 
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Figure 2-5. Time-Series plot options page (From: www.oscar.noaa.gov) . 
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Data Display & Download 


Data Download 



Time^Longitude Section for Ocean Surface Current 


Non-Java Version 


Plot Type: Zonal and Meridional Currents 
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Figure 2-6. Time-Longitude section plot options page (From: www.oscar.noaa.gov) . 

D. DATA DOWNLOAD AND DISPLAY 

Data can be downloaded from the website. Geographical limits and other 
parameters can be applied through the Data Download interface (Figure 2-7). The 
download interface gives the user a variety of options to download data per one’s unique 
needs. For download only, unfiltered data is available and the user can select 5-day mean, 
monthly mean or long-term mean data sets. In terms of variables, either (U, V) mean or 
anomaly data can be downloaded. All data sets can also be downloaded from the dapper 
server (Sirott et ah, 2004) at http://dapper.pmel.noaa.gov/dapper/oscar . All the 
downloadable files are in the netCDF format and can be further processed with Matlab. 


13 





















Data Display & Download 



Data Download 

Lat-Lon Map 

Lat-Time Section 

Time Series Plot 

Time-Lon Section 


Data Download for Ocean Surface Current NonOava Version 


Data Type: @ 5-day Mean © Monthly Mean © Long Term Mean (1993-2006) 

Filter Type: @ unfilered 
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Figure 2-7. Data Download options page (From: www.oscar.noaa.gov) . 

One of the handiest features of OSCAR website is that a user can directly get 
plots for the region of interest. Online surface currents’ plots can be generated in the 
shape of a Lat-Lon Map, Lat-Time Section, Time-Lon Section, and Time-Series Plot. For 
online plots, a variety of options are available to generate the desired surface currents’ 
plot. Selectable options are “Data Type” in which one can select 5-day interval data, 
monthly mean, mean over selected range and long- term or seasonal mean (Figure 2-4). 
“Plot Type” describes the shape of the final plot that can be a vector mean or a vector 
anomaly plot. Plots can also be generated solely for the meridional or zonal current 
components (Figures 2-9 and 2-10). Users can also set the time and geographical limits of 
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the final plot. For online display, both filtered and unfiltered data are available. For 
veetor overlay anomaly red eolor has already been set in the system for eastward anomaly 
(Figure 2-8) and blue for westward anomaly. Sealing ean also be applied to the plots. 

Time-Series plot (Figure 2-11) ean be generated for any desired range of time 
starting from 1992 either with filtered or unfiltered data. Also Lat-Time plot (Figure 2- 
12) and Time-Lon plot (Figure 2-13) can be generated online for the desired selection of 
Latitude and Time or Time and Longitude. 
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Figure 2-8. (Top) Vector overlay of surface currents; (Bottom) Vector overlay of 
surface currents anomaly. Large eastward anomalies can be seen in eastern 
equatorial pacific (From: www.oscar.noaa.gov) . 
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6-Day Interval Ocean Surface Currents (meter/sec) 

Centered on May 3 2008 
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Figure 2-9. Latitude-Longitude map of Meridional and Zonal Current Mean centered 
on May 2, 2008 (From: www.oscar.noaa.gov) . 
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6-Day Interval Ocean Surface Currents (meter/seo) 

Centered on May 2 2008 
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Figure 2-10. Latitude-Longitude map of Meridional and Zonal current anomaly 
centered on May 2, 2008 (From: www.oscar.noaa.gov) . 
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Time Series Averaged Over Region: 30.CIE-380.5E, 60.6S—60.5N 


Figure 2-11. 


Figure 2-12. 
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Time-Series plot (From: www.oscar.noaa.gov) . 


5—Day Interval Ocean Surface Currents at 140'’¥ 


eo'»N 

5a*N 

4a'»N 

3Q'»N 

30'»N 

ia*N 

0 " 

lO^S 

S0"S 

30*S 

50*3 

60"S 



I o.a 

■ 0.4 

■ 0.0 

■ -0.4 

t - 0.8 
-l.S 


60'»N 

5a*N 

4a'»N 

3Q'»N 

SO-N 

10*N 

0“ 

lO^S 

20'^S 

30*S 

40’S 

50’S 

60’S 



0.1 

0.0 


-o.i 


1992 1995 1998 


2001 2004 2007 


NE£DIS/l4CM 


Noy 7 3oaa 


Latitude-Time section of Zonal and Meridional currents (From: 
www.oscar.noaa.gov) . 
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Figure 2-13. 


Time-Longitude seetion of Zonal and Meridional eurrents (From: 
www.oscar.noaa.gov) . 
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III. OPTIMAL SPECTRAL DECOMPOSITION (OSD) METHOD 


A. INTRODUCTION 

A major weakness of the OSCAR dataset is its inability to reveal the coastal 
boundary currents, such as the Gulf Stream and the Kuroshio. Besides, the data is quite 
noisy (Figure 3-1). 


Of|okwlD«tt2007 Jan A 



Longliidft 


Figure 3-1. Global surface velocity vectors on January 2-6, 2007 (From: OSCAR 
data) 

To improve the deficiencies of the OSCAR data system, this study uses the 
Optimal Spectral Decomposition (OSD) method. Any field (temperature, salinity, or 
velocity) can be decomposed into generalized Fourier series using the OSD method. The 
three-dimensional field is then represented by linear combination of the products of basis 
functions (or the so-called modes) and corresponding Fourier coefficients. If a 
rectangular closed ocean basin is considered, the basis functions are sinusoidal functions. 
If a realistic ocean basin is considered, the basis functions are the eigen-values of the 
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three-dimensional Laplace operator with real topography. The Fourier coefficients are 
determined from observational data through solving a set of linear algebraic equations. 
One major benefit of using the OSD method is that the boundary conditions for the ocean 
variables (temperature, salinity, and velocity) are always satisfied (Chu et ah, 2003a, b; 
2004, 2007). 

B. SPECTRAL DECOMPOSITION 

Following Chu et ah, 2003a, and Chu et ah, 2005, a velocity field U(x^,0 at any 
point of the area of interest fl can be written in the form of parameter-weighted sums of 
the harmonic Z^(x^) and basis 'F^(xq) functions by 

U(Xo,0 = X^*<^^)[*^^'^^^(*o)] + Z^-t(0[kxV'F^(Xo)] (3-1) 

s=l k=l 

The coefficients b^(t) and a^(0 in Eq. (3-1) are functions of time; k is the unit 
vector orthogonal to the surface z=const and pointing upward; and V = (V^,V2) is the 
horizontal gradient operator. A set of harmonic functions, , is introduced to account for 

the multi-connection feature of S (i.e., existence of islands) and/or disconnected rigid 
boundaries (if applicable). 

The harmonic and basis functions used in (3-1) satisfy the following equations 

Z1Z,=0, (3-2) 

where A is the horizontal Laplacian operator and are eigen-values of corresponding 

spectral problem. Boundary conditions (to be specified in the next Section), in general, 
assume imposing of no-flux through a rigid boundary and normal derivatives along an 
open boundary. 

For the given harmonic, Z; (x^), and basis ,'F^(xq) , functions, and a series of 
original OSCAR velocities, , Eq. (3-1) is used to construct an appropriate variation 
problem and estimate unknown spectral coefficients a = [aj,...,a^]', , and 

the optimal truncation number K. Here, is the transpose operator. When the spectral 
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coefficients, a = [ap...,a^]', b = , c = [cj,...,c^ and the optimal truncation 

numbers, K and M , have been estimated, Eq. (3-1) is used again to project the velocity 
on the OSCAR grid. 


C. OPTIMAL MODE TRUNCATION 

The optimal mode truncation number {Mopt) is defined as the critical mode 
number with the set of spectral coefficients {ak} least sensitive to observational data 
sampling and noise. For sample size of P and mode truncation of M, the spectral 
coefficients [ak} are estimated by the least square difference between observed and 
calculated values, 

(J) 

where the symbol represents the estimated value at (x, t). For homogeneously 
sampled data, with low noise and without systematic error, the empirical cost function 
should tend to 0 monotonically as M increases to infinity. The set of the spectral 

coefficients [uk] depends on the mode truncation M. Optimal estimation of [ak] is 
equivalent to the determination of Mopt [Ivanov et ah, 2001a and Chu et ah, 2003 a, b]. 


1 

■^emp ~ ^ 


^min 



is used to determine the optimal mode truncation Mopt- Here, r is the probability of 
J-J ^0 asMincreases. 


For the global oceans, Mopt is 30-50 for the basin-scale (-300 km) variability and 
150 for the mesoscale (-20 km) variability [Chu et ah, 2005a]. For sparse and noisy data, 
it is difficult to get reliable and stable estimates of all the necessary spectral coefficients, 
but the first few spectral coefficients are reliable and stable. Figure 3-2 shows the global 
velocity vectors on January 2-6, 2007 after using the OSD method on the OSCAR data. It 
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shows that the global circulations are much more reasonable with evident boundary 
currents, such as the Gulf Stream and the Kuroshio. 

O8D«IIID0ihlMtB2(n7 JOD 4 



Longlfcjds 

Figure 3-2. Reconstructed global surface velocity vectors on January 2-6, 2007 (From 
OSCAR data using the OSD method). 
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IV. SEASONAL VARIABILITY 


A. DATA PROCESSING 

The monthly mean surface currents were reconstructed from the OSD-reprocessed 
OSCAR data. This data was derived from 1993 - 2006 satellite altimetry. This study will 
first look at the whole Indian Ocean during the two monsoon regimes; then it will 
discuss, in detail, the circulation patterns of the three most dynamic regions of the Indian 
Ocean (Somali Current, North Arabian Sea, and Bay of Bengal). 

B. SURFACE CIRCULATIONS 

In Chapter I, the study showed visuals of the conventional seasonal wind flow in 
the Indian Ocean region, which is the main contributing factor in driving the ocean 
surface circulations. In this chapter, this study will analyze the seasonal variability in 
surface circulation patterns for the Indian Ocean using the OSD-reprocessed satellite 
altimeter data. For the whole Indian Ocean, this study analyzes the two typical monsoon 
regimes and performs side-by-side comparisons with results from the most recent 
studies. 


C. WHOLE INDIAN OCEAN 

I. Southern Hemisphere 

Figures 4-2 and 4-4 show a broad zonal eastward flow of the South Equatorial 
Current (SEC). The southeasterly Trades drive this (Eigure 4-1). This broad zonal flow 
supplies the western boundary currents along the eastern coast of Madagascar (Eigure 4- 
2). The same flow pattern has also been observed by Schott and McCreary, (2001). In 
about the center of Madagascar, the SEC splits into northward and southward flowing 
branches. The northern branch, which is known as the Northeast Madagascar Current 
(NEMC), supplies the East African Coast Current (EACC). The southward branch flow 
as the Southeast Madagascar Current (SEMC). Erom the results produced by the 
altimetry data, it was observed that ocean currents south of about 15° S do not show 
much seasonal variability. Therefore, results from altimeter data, Eigures 4-2 and 4-4 
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show almost the same branches in SEC regimes in both seasons. The same has been 
observed by others e.g., Swallow et al. (1988), and represented in an earlier study by 
Schott and McCreary (2001) in Figures 4-3 and 4-5. 


2. Summer Monsoon 


Figure 4-2 shows the mean surface circulation of the Indian Ocean during July 
derived from the satellite-based altimeter data. During the summer monsoon as described 
in the previous studies, the SEC and EACC supply the northward flowing Somali current. 
Figure 4-2 clearly shows that part of Somali current is recirculating to form the “Southern 
Gyre” near the south of the equator. The other part is moving in the northward direction. 
A gyre, known as the “Great Whirl,” is formed” just north after crossing the equator. The 
same has also been observed by Schott and McCreary (2001). Figure 4-3 also shows the 
presence of a small eddy, the “Socotra Eddy,” which is not evident in the current field 



Figure 4-1. Monsoon wind stress from the NCEP climatology for a) January; b) April; 
c) July; d) October (From: Schott & McCreary, 2001). 
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During summer, the Southwest Monsoon Current (SMC) flows eastward (Figure. 
4-2). This flow is also confirmed by Schott and McCreary (2001) (Figure 4-3). Figure 4-2 
shows that SMC is supplied by the outflow of Somali current and that some of its supply 
seems to be coming from the southward-flowing West Indian Coastal Current (WICC). 
Along the eastern coast of India, we can see a weak flow of the East Indian Coastal 
Current (EICC) which flows northeastward all along the coast (Eigures 4-2 and 4-3). The 
EICC gets most of its supply from the Southwest Monsoon Current (SMC) to the east of 
Sri Lanka. Along the eastern coast of Madagascar, South Equatorial current (SEC) splits 
up into two branches: one moves northward as the Northeast Madagascar Current 
(NEMC); the other moves southward as the Southeast Madagascar Current (SEMC). 
Presence of an eddy can also be seen towards the northwestern edge of Madagascar. This 
has also confirmed by the previous studies. In case of the Leeuwin Current, the general 
southward flowing direction cannot be seen in the results of this study. This is a weakness 
of the satellite altimetry data. 


Jul 



20“^ 40*^ eo'fe SO't 10Ot 120^ 


Longitude 

Eigure 4-2. Mean Surface Circulation pattern of Indian Ocean during July, derived 
from satellite-based altimeter data. 
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Figure 4-3. A schematic representation of identified current branches during the 
southwest monsoon, including some chokepoint transport numbers 
(Sv=10^6m^3s^-1). Current branches indicated (see also Figure 9) are the 
South Equatorial Current (SEC), South Equatorial Countercurrent 
(SECC), Northeast and Southeast Madagascar Current (NEMC and 
SEMC), East African Coast Current (EACC), Somali Current (SC), 
Southern Gyre (SG), and Great Whirl (GW); and associated upwelling 
wedges: Socotra Eddy (SE), Ras al Hadd Jet (RHJ), and upwelling wedges 
off Oman, West Indian Coast Current (WICC), Laccadive High and Low 
(EH and EL), East Indian Coast Current (EICC), Southwest and Northeast 
Monsoon Current (SMC and NMC), South Java Current (JC), and 
Leeuwin Current (EC). See text for details. (Erom: Schott & McCreary, 
2001 ). 


3. Winter Monsoon 


Eigure 4.4 shows the mean surface circulation of the Indian Ocean during 
December, derived from the satellite based altimeter data. It can be seen that the 
northward flowing East African Coastal Current (EACC) and the southward flowing 
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Somali Current meet in a confluence zone at about 1-4° S. Then both supply the eastward 
flowing South Equatorial Countercurrent (SECC) (Eigure 4-4). This pattern is the same in 
Schott and McCreary (2001) (Eigure 4-5). On the eastern end of the SECC, the South 
Java Current (Eigure 4-4) flows southeastward. The same circulation pattern of SECC 
can also be seen in Eigure 4-5. 


Dec 
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Longitude 

Eigure 4-4. Mean Surface Circulation pattern of Indian Ocean during December, 
derived from satellite-based altimeter data. 

During this time of the year, the southerly flowing Somali current is supplied 
from the westward flow in the Arabian Sea. This inflow into the Somali current seems to 
be coming in several branches. A lower latitude supply is originating from the westward 
flowing Northeast Monsoon Current (NMC) south of Sri Eanka. This also supplies the 
WICC after circulating around the Eaccadive High. This circulation feature is not 
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prominent in the results generated for December, but it can be seen as further in winter 
(Figure 4-6). This feature was also explained by Bruce et al. (1994). The second branch 
that seems to be inputting the Somali current comes from the westward flow south of 
Socotra. This flow maybe getting input from WICC as well (Figure 4-4). 



Figure 4-5. As in Figure 4-3, but for Northeast Monsoon (From: Schott & McCreary, 
2001 ) 

After analyzing the circulation pattern in the Bay of Bengal, this study observed 
that the EICC reverses direction twice a year: (Figures 4-2 and 4-7a) when it is moving 
northeastward and, after that, (Figure 4-7b) when it is moving in southwestward. In the 
interior of the Bay, a basin-wide clockwise gyre (Figure 4-7a) was observed from 
February to May. This gyre disappears during the summer monsoon, when the boundary 
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Figure 4-6. Mean Surface Circulation pattern of Indian Ocean during January, derived 
from satellite-based altimeter data after OSD. 

4. Equatorial Regime 

A phenomenon, which is unique to the Indian Ocean due to the semiannual 
eastward winds along the equator (Figure 4-lb,d), is the flow of strong eastward surface 
jets (Wyrtki, 1973). These jets appear during the transition seasons between the 
monsoons, i.e., April to June and October to December. These jets carry warm upper- 
layer waters eastward. They lower sea levels and decrease mixed-layer thickness in the 
west, but increase them in the east (e.g., Rao et ah, 1989; Figure 4-8b,d), thereby causing 
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current splits at the 10°N confluence as observed by Schott and McCreary (2001). During 
the winter season, there seems to be a clockwise circulating gyre in the interior of the 
Bay. The EICC attached to it seems to have a southward flow. 
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an eastern temperature maximum in the 100-m map (Figure 4-8d,f). The strong eastward 
flow of these jets can be seen in Figures 4-7a and b along the equator. These eastward¬ 
flowing, upper-layer, current anomalies propagate pole-ward along the eastern boundary 
of the basin (Schott and McCreary, 2001). 
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Figure 4-7. Mean Surface Circulation pattern of Indian Ocean during a) April, 
October; derived from satellite altimeter data. 
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Figure 4-8. Climatological temperature at a-b) The surface during January and July, 
c-f) 100 m depth during January, May, July, and November, areas of 
relatively high temperature are shaded (From: Rao and Sivakumar, 2000). 

D. SOMALI CURRENT 

Traditionally, analyzing the annual cycle of the Somali Current system with its 
reversal of an entire current system (Figure 4-9) has remained an area of strong interest to 
the Indian Ocean oceanographic community. This section briefly analyzes the seasonal 
variability of the Somali Current and also compares this study’s results with previous 
studies on this topic. 


31 


















Longitude Longitude 


Figure 4-9. Reversal of Somali Current derived from mean altimeter data. 

1. Seasonal Development 

To see the seasonal development of the Somali Current, data for the specific area 
was extracted and vector plots were made for each month. The following paragraphs will 
show how the Somali current develops seasonally and compare this study’s findings side- 
by-side with published studies. Figure 4-13 shows the seasonal development of the 
Somali Current system originally described by Schott et ah, (1990) and was updated by 
the same author in 2001. 

March-May: From satellite altimeter-derived data, this study observed that, 
before the onset of the monsoon, the southern Somali Current is an extension of the East 
African Coast Current (EACC) that flows northward. Just about south of the equator it 
turns offshore. The data also shows evidence of a gyre, which is most likely the Great 


32 






































Whirl. Previous studies do not show the presence of such a gyre during this time of the 
season. At this time the season width scale of the Somali current is of the order of 50-100 
km (Schott and McCreary, 2001). 


Apr Jul 



Longitude Longitude 


Figure 4-10. Mean surface Somali currents for a) April, b) July, derived from satellite- 
based altimeter data. 

June-July: As the monsoon season starts in June, the “Great Whirl” is seen on 
lower latitude than in April (Figure 4-10b). Schott and Quadfasel (1982) made 
observations from moorings in the northern Somali Current regime. During their 
observational year, there was an unexpected start of the monsoon and, after this start, 
they reported distinct westward-propagating signals. They understood them as the first¬ 
mode of Rossby waves. Their conclusion was that the start of the Great Whirl was a 
response to the strong anticyclonic wind-stress curl offshore of the Somali coast by these 
long Rossby waves. They reflected into short Rossby waves at the boundary and 
accumulated energy there. At this time of the year, the Somali Current is strongest and 
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tends to flow northward along the coast. The observation made by Schott et ah, (1990) 
also explains the creation of a gyre at about 4°N, where part turns eastward and part flows 
back across the equator in a circulation pattern which is termed as the ‘Southern Gyre.” 
This also shows up in this thesis study’s results, but at lower latitude at about 4°S (Figure 
4-10b). Why it shows up at a lower latitude needs more study which has not been covered 
in this study. 

August-September: In the after part of the summer monsoon, the Great Whirl 
seems to be maintaining its position (Figure 4-11 a). Transports in this late summer 
monsoon phase can exceed 70 Sv (Fischer et ah, 1996; Schott et ah, 1997). During this 
time, the Somali current flows generally northward all along the coast. 

October-November: At the end of the southwest monsoon season, the Somali 
Current turns offshore at 2-3°S, while the Great Whirl has now shifted to a higher latitude 
(Figure 4-1 lb). During this time of the year, the flow is generally northward but its 
strength has reduced definitely. The evidence of Great Whirl even present till the end of 
the year is also explained by Bruce et al. (1981). 



Figure 4-11. Mean Somali current circulation pattern during a) August, b) October, 
derived from satellite altimeter data. 
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Winter Monsoon 


In the winter northeast monsoon, the winds blow southward away from the Indian 
subcontinent (Figure 4-la). Due to this reversal of the wind pattern, the Somali Current 
also reverses its flow towards the south. During winter, the southward flowing cross- 
equatorial Somali current meets the northward flowing EACC in a confluence at about 2- 
4°S. Then, both flows input the eastward flowing South Equatorial Countercurrent 
(Eigure 4.12). This flow pattern is also confirmed by Du'eing & Schott (1978) and 
Swallow et al. (1991). At the equator, the southward Somali Current is quite shallow, 
carrying 5 Sv in the upper 150 m. This is because there is a northward undercurrent at 
this time (Schott and McCreary, 2001). 


Jan 
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Eigure 4-12. Mean Somali current circulation pattern during winter monsoon, derived 
from satellite-based altimeter data. 


Even at this time of the year, the Great Whirl is seen between 0-5°N. Numerous 


observations had been made in the northern Somali Current for the winter monsoon. 
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Bruce et al. (1981) reported that the Great Whirl can actually sustain summer monsoon 
conditions below the developing winter monsoon surface circulation for some time into 
the winter monsoon. 

SOUAU CURRENT FLOW PATTERNS 



ap|»rl^rK(-t-1«0<i4 

indHonenls (- IW - MOm) 


Figure 4-13. Schematic diagram of Somali current upper layer flow pattern over the 
course of the year (From: Schott & McCreary, 2001). 


SOMALI CURRENT PLOW PATTERNS 



Figure 4-14. Schematic Diagram of Somali current upper layer pattern over the course 
of the year derived from satellite-based altimeter data. 
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3. Comparison with Earlier Studies 

Comparing this study’s results (Figure 4-14) from satellite altimeter data with 
Schott and McCreary (2001)’s schematic diagram (Figure 4-13) from ship drift data, it 
can be observed that the general circulation pattern of the Somali Current is almost the 
same. Few features are exceptions, such as the position of the southern gyre and the true 
depiction of the currents near the land boundaries. Thus, with the altimeter’s data, this 
study clearly observed the mesoscale features of a certain part of the ocean. However, it 
is difficult to distinguish between the small-scale features and the features along the 
coast. 

E. NORTH ARABIAN SEA 

In the North Arabian Sea, broad eastward flow during the summer monsoon 
across the whole latitude range of the Arabian Sea can be seen (Figure 4-15c). The exact 
condition is also observed from the ship drift currents (Figure 4-16a). In this study’s 
results, we cannot see the presence of the Laccadive Low off the southwest tip of India 
(Figure 4-15c) during the summer Monsoon season. Schott and McCreary (2001) 
observed this in previous studies. As we move further towards the end of summer 
monsoon in October (Figure 4-I5d), one can clearly observe the circulation of the 
Laccadive Low along the southwestern tip of India moving in an anti-clockwise 
direction. At this time of the year, the SECC has started flowing eastward almost along 
the equator. The southwest monsoon current flows eastward along the southern tip of 
India. 

As part of the Joint Global Ocean Flux Study (JGOFS), for the years 1994-1996, 
the northern Arabian Sea has been comprehensively studied. In the JGOFS study during 
the summer monsoon, it was observed that a northeastward coastal jet developed along 
the south coast of Oman in early May which persisted throughout the summer. This jet 
has been called Ras A1 Hadd Jet by Bo"hm et al. (1999). This jet flows past Ras al Hadd 
at the southeast comer of Oman (Figure 4-3) and meanders into the Gulf of Oman. As 
explained by Flagg and Kim (1998) and Bo"hm et al. (1999) this jet is linked to an 
anticyclonic gyre just south of its offshore turning point and a cyclonic gyre to its north. 
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These small-scale gyres cannot be seen in this study, but the flow pattern of the jet can be 
seen along the eastern tip of the Arabian Peninsula (Figure 4-15c). 



Longitude Longitude 

Figure 4-15. Mean surface circulation pattern of the North Arabian Sea all through the 
year, derived from satellite-based altimeter data, showing a) January, b) 
April, c) July, d) October. 
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Figure 4-16. Surface current vectors from drifter climatology; a) July - September, and 
b) December - February (From: Schott and McCreary, 2001). 

When observing the winter circulation pattern in the North Arabian Sea (Figure 4- 
15a), a basin wide clockwise gyre and, along the northern coast, weak easterly currents 
were seen. As we move down towards the south, westerly currents appeared that can also 
be seen in the ship drifts in the northern Arabian Sea (Figure 4-16b). In January and 
April, the western-boundary current of the Gulf of Oman north of Ras al Hadd is 
comparatively weak which is also observed by Flagg & Kim (1998). As explained by 
Schott and McCreary (2001) (Figure 4.3), the presence of the Laccadive High (LH) can 
be seen along the southwestern tip of India moving in the clockwise direction. The SECC 
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is flowing eastward and the northeast monsoon current is flowing in the westerly 
direction. This has also been observed by Schott & McCreary, 2001 (Figure 4-5). In 
general during the winter monsoon, the mean flow is westward. Shi et ah, (2000) also 
explains this with the addition that this westward flow is accompanied by down-welling 
along the coast and an offshore-directed pressure gradient. Moving further offshore, a 
mesoscale eddy along the equator at about 71° E (Figure 4.15a) moving in anti-clockwise 
direction is seen. This mesosacle feature needs more study because studies carried out 
during 2001 also do not explain anything about it. 

As the winter northeast monsoon dies down in April (Figure 4-15b), the general 
circulation pattern is still the same in the northern Arabian Sea with some weak westward 
flow. Along the coast of Pakistan we can see weak eastward current. During this time of 
the year the Somali Current has started moving in northeasterly direction along the coast. 
At about equator it comes in confluence with the East African Coast Current (EACC) 
which is supplying the South Equatorial Counter Current (SECC). The Great Whirl can 
also be clearly seen during this time of year. 

F. BAY OF BENGAL 

The surface circulations in the western part of the Bay of Bengal are strongly 
influenced by the seasonal reversal of the wind stress forcing. The climatological 
situation in the western part of the bay is similar to that in the western Arabian Sea. For 
EICC, the local winds have weaker effects because remote forcing is relatively more 
important for the EICC than they are for the Somali Current. As explained by Schott and 
McCreary (2001), the crossing of remotely forced coastally trapped-waves, and Rossby 
waves in the interior of the bay, are major contributors to the EICC variability. One of the 
other major factors, which prevents the accurate determination of only geostrophic 
currents in the bay, is the large fresh-water flux from both rainfall and river discharges 
(Schott & McCreary, 2001). 

The seasonal cycle of the surface circulation in the Bay of Bengal from satellite 
altimetry is presented in Figure 4-17. We observe that as for the Arabian Sea, circulation 
patterns in the interior of the Bay of Bengal contain many mesoscale features, but there is 
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no evidence of any seasonally reoccurring coastal gyres like the Southern Gyre and the 
Great Whirl. As the summer monsoon begins, the northward surface currents in the EICC 
off India do not strengthen. However, there are weak northward currents along the 
eastern coast of India (Figure 4-I7c). In the summer, there are hardly any mesoscale 
features inside the Bay of Bengal. Rather, the general flow inside the bay is eastward at 
this time of year. The same pattern can also be seen by results from ship drift data shown 
in Figure 4-18a. This eastward flow is weaker along the coast and strengthens as it flows 
further south inside the bay. Its maximum strength is located just south of the bay. That is 
the southwest monsoon current. In October (Figure 4-I7d) when the summer monsoon 
has almost diminished, an entirely reversed circulation pattern inside the bay is present. 
Now EICC seems to be moving southward along the east coast of India and has more 
strength as compared to it in summer. A mesoscale feature can also be seen inside the 
Bay as a large-scale eddy appears in the southwestern boundary of the Bay along EICC. 
In the eastern part of the Bay, the weak currents are flowing in a northerly direction. 
Inside the bay, the currents are moving mostly counterclockwise. 

During winter monsoon in January (Figure 4-17a), circulation inside the Bay is 
completely transformed into a mesoscale eddy moving in a clockwise direction. The 
strength of the currents inside the Bay has also increased as compared to summer. In the 
southern part of EICC, confluence occurs at about 17° N. Here another smaller eddy is 
moving in an opposite direction. Moving further south the general flow is towards the 
west. In an analysis of sections observed during December 1991, between the northern 
end of the Bay and 10°N, Shetye et al. (1996) found the geostrophic transport relative to 
1000 dbar increased southward from 2 Sv to about 8 Sv. 

The winter monsoon dies down in April. The clockwise circulation pattern inside 
the bay still persists while the center of circulation now has shifted towards western side 
of the bay (Figure 4-17b). During this time of year, EICC seems to be flowing towards 
north along the east coast of India. In the southern part of the bay, the general circulation 
is towards the west. This circulation pattern is also confirmed by Shetye et al., (1993) 
who also estimated a transport of about 10 Sv for the northward EICC at this time and 


41 



Latitude ^ Latitude 


concluded that it was the western-boundary current associated with a wind-driven 
anticyclonic gyre in the interior of the bay. 
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Figure 4-17. Mean surface circulation pattern of the Bay of Bengal all through the year, 
derived from satellite-based altimeter data showing a) January, b) April, c) 
July, d) October. 
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Figure 4-18. Surface currents in Bay of Bengal from the ship drift data; a) July, b) 
November (From: Eigenheer & Quadfasel, 2000). 

G. MONSOON CURRENTS 

In this subsection, the focus is on the variability solely in the monsoon currents in 
the northern Arabian Sea. The region between Sri Lanka and the equator is a bottleneck 
for the zonal currents in the northern Indian Ocean. They provide the major pathway for 
the exchange of water between the Arabian Sea and Bay of Bengal. As known from ship 
drift observations (Cutler and Swallow, 1984), this current reverses with the monsoons. It 
flows eastward as SMC from June to September (Figures 4-15c and 4-16a) and westward, 
as NMC, from November to February (Figures 4-15a and 4-16b). 
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Figure 4-19. Monsoon current transport south of India - Sri Lanka; a) Decomposition 
into depth bins using two different extrapolation schemes from top 
instrument level towards the surface: extrapolation of moored current 
meter shear from upper most levels (solid) or using ship drifts as surface 
values (dashed). The heavy curve on the left hand side shows annual 
mean, b) transport time-series, based on shear extrapolation; dots mark 
shipboard ADCP sections for comparison (From: Schott et ah, 1994). 

1. Southwest Monsoon Current (SMC) 

Figure 4-2 shows the flow pattern of the southwest monsoon current at the 
southern tip of India. During the summer, this current moves eastward with greater 
strength than other current regimes near it (Figures 4-2 and 4-3). The same flow pattern 
of SMC was also observed by Schott et al. (1994). The SMC decayed so quickly with 
depth that its transport had to be determined either by extrapolating ADCP and current 
meter observations to the surface or by using seasonal mean ship drift currents as surface 
values (Figure 4-19a). Based on shear extrapolation, the eastward SMC transport between 
4-6°N was estimated to be 8.4 Sv, with 70% in the top 100 m. From the current and 
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transport time-series, they also saw the strong intraseasonal variability (Figure 4-19b). 
The source of water for the SMC eventually must be in the western Arabian Sea. After 
passing Sri Lanka, much of the SMC turns northward into the Bay of Bengal; whereas 
some of it turns southward (Figure 4-2). Finally, this current inputs to the current systems 
along the western coast of Malaysian islands. 

2. Northeast Monsoon Current (NMC) 

The NMC flows westward south of India (Figures 4-4 and 4-5). During 1991, its 
transport was observed to be 11 Sv (Schott et ah, 1994). This transport is larger than the 
estimated SMC transport and results in a net annual-mean transport of 2.5 Sv to the west 
for 1991 (Figure 4-20a). In this study’s results, the NMC flow persists through March. 
Hacker et al. (1998) saw the same. The NMC partly carries fresher water from the Bay of 
Bengal the properties of which are detectable in the WICC (Stramma et ah, 1996). It 
seems that most of the source of NMC is from the recirculating water brought eastward 
north of the equator. This finding is the same as concluded by Hacker et al. (1998). 
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V. INDIAN OCEAN DIPOLE AND EASTWARD 
EQUATORIAL JET 


A. THE INDIAN OCEAN DIPOLE (lOD) 

In 1999, the Indian Ocean Dipole (lOD) was discovered by Professor Yamagata, 
Dr Saji, and associates of the Climate Variations Program of Frontier Research System 
for Global Change (Vinayachandran et ah, 1999). They described it as an inter-annually 
occurring climate mode in the tropical parts of the Indian Ocean. The occurrence of 
Dipole Mode events can be traced back through the coral record to the mid-Holocene 
period which dates about 11,000 years before present to present. As typical to other 
climate oscillations, the lOD climate mode has a ‘positive’ phase and a ‘negative’ phase 
(Figures 5-1 and 5-2). 


Positive Dipole Mode 



Figure 5-1. Schematic diagram of SST anomalies (red shading, warming; blue, 
cooling) during a positive lOD event. White patches indicate increased 
convective activity. Arrows indicate wind direction. (From: 
A.Suryachandra Rao. Frontier Research Center for Global Change). 
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During a positive lOD event (Figure 5-1), the sea surface temperature (SST) rises 
in the western equatorial Indian Ocean, which is along the eastern coast of Africa, and 
from the northern half of Madagascar to the northern edge of Somalia. The SST drops in 
the southeastern part of the Indian Ocean, which is along the northern coast of Australia, 
the eastern coast of Japan, and throughout Indonesia. In Figures 5-1 and 5-2, the red color 
denotes rise in SST and the blue color corresponds to fall in SST. Associated with these 
changes, the normal convection situated over the eastern Indian Ocean warm pool shifts 
to the west and brings heavy rainfall over eastern Africa and severe droughts/forest fires 
over the Indonesian region. 


Negative Dipole Mode 



Figure 5-2. Schematic diagram of SST anomalies (red shading, warming; blue, 
cooling) during a negative lOD event. White patches indicate increased 
convective activity. Arrows indicate wind direction. (From: 
A.Suryachandra Rao., Frontier Research Center for Global Change). 

Prof. Yamagata, Dr. Saji, and other researchers of the climate variations research 
program of Frontier Research Center for Global Change (FRCGC) created the name of 
lOD to represent the zonal dipole structure of the various coupled ocean-atmosphere 
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parameters, such as SST, OLR, and Sea Surface Height anomalies. Furthermore, 
convective patterns increase in the northern half of Africa, India, and off the eastern coast 
of Africa. Inverse conditions exist during a negative lOD event (Figure 5-2). 

B. DIPOLE MODE INDEX 

The lOD events are represented by the Dipole Mode Index (DMI), which is 
defined as the zonal difference of sea surface temperature anomaly between the tropical 
western Indian Ocean (50°E-70°E, I0°S-I0‘’N) and the southeastern Indian Ocean (90‘’E- 
1I0°E, lO^S-equator) (http://ioc3.unesco.org/oopc/state of the ocean/sur/ind/dmi.php) . 
The time series of the monthly mean DMI (Eigure 5-3) clearly reveals the evidence of a 
large positive anomaly during 1994, 1997, and 2006. The most prominent amongst them 
is during 1997-1998. Although negative anomalies are not as strong as positive ones, they 
can be observed during years 1996, 1998, 2001, and 2005. The most prominent amongst 
them is 1998. Eor the subsequent analysis, the positive lOD events of 1997-1998 and the 
recent one during 2006 are considered. To see our results in context to the negative lOD 
events, the concentration will be on the events during 1998-1999 and 2001. 



Eigure 5-3. Dipole Mode Index for the period 1992 - 2006. 
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C. THE EASTWARD EQUATORIAL JET (OR WYRTKI JET) 


One of the most prominent features in the Indian Ocean Equatorial region is the 
eastward flowing Equatorial jets that usually can be observed during the transition period 
between monsoons (Wyrtki, 1973). Several modeling studies were carried out to study 
the connection between the eastward equatorial jet and the lOD events. The results 
showed that the eastward equatorial jet weakens in a positive Dipole Mode event and 
strengthens in a negative Dipole Mode event (Vinayachandran et ah, 2001). This feature 
has also been observed by Reppin et al. (1999) while making observations in the Indian 
Equatorial region during 1994’s Positive Dipole event. This shows that, in context to the 
surface circulation in the Indian Ocean, the behavior of these jets tend to vary more than 
any other large-scale feature during any Dipole Mode event. 

Satellite-based observational surface currents data was used to investigate the 
long-term connection between the eastward equatorial jet and the lOD events. Therefore, 
for the subsequent analysis only the behavior of these jets in response to different Dipole 
Mode events will be viewed. These eastward flowing jets can be seen along the equator 
in Eigures 4-7a and b. 

D. COMPLEX EFO ANALYSIS OF THE CURRENT DATA ALONG THE 

EQUATOR AND COMPARISON WITH DMI 

To see the variance in the spatial patterns along the equator from 45° E to 95° E, 
data was extracted for the area of interest. Then, a calculation of the 13-year arithmetic 
mean of the data set by using the basic statistical relation was performed, 

Z n 

_ i=\ 

n n 

Complex Empirical Orthogonal functions were used to calculate the amplitude of 
the first four modes. The variance was calculated by computer-based program. Therefore, 
there were a finite number of modes that could be derived. In our case it was 50. After 
calculating, it was found that the first four modes cover the 91% of the variance as 
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illustrated in Table 5-1. In this case, to see the variance in surface circulation pattern, 
91% variance could have revealed the bigger anomalies in the spatial patterns. Therefore, 
only the first four modes were taken into account for the subsequent analysis. 


EOF Mode 

Relative Variance 

Cumulative Relative 
Variance 

Mode 

66.95 

66.95 

2"^ Mode 

13.27 

80.22 

3'^" Mode 

6.98 

87.19 

4“' Mode 

4.23 

91.42 


Table 5-1. Percentage of Variance of first four EOF Modes. 

The Amplitude of the Variance for all four modes was analyzed to see if any 
relation exists with the general trend in DMI. 


Amplitude Mode 1 



Figure 5-4. Amplitude of EOF Mode-1. 
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Amplitude Mode 2 




Figure 5-5. Amplitude of EOF a) Mode-2, b) Mode-3. 
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Ampltude Mcde 4 



Figure 5-6. Amplitude of EOF Mode-4. 

Figure 5-4 represents the amplitude of the first EOF mode. While comparing the 
amplitude of Mode-1 with the general trend in DMI, it was seen that the biggest variance 
is analogous to the trend in DMI during the year 1997-1998, the year of strongest positive 
Dipole Mode event in the recent history. Although not very conspicuous, small peaks of 
variance can be seen in the spatial patterns of currents along the equator during the years 
1994 and 2006. These are again comparable with DMI and are the years of positive 
Dipole Mode event. 

When amplitudes of variance from mode-2 to mode-4 were compared, it was 
clearly seen that the biggest anomaly during 1997-1998 was evident in all three modes. In 
mode-2, the evidence of a negative mode event that occurred during year 2005-2006 was 
also viewed. Therefore, generally speaking, the spatial anomaly of the currents along the 
equator to the occurrence of a Dipole Mode event can be directly related. 
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E. EQUATORIAL JET DURING POSITIVE DIPOLE MODE EVENTS 

To analyze the behavior of the equatorial jets during known positive Dipole Mode 
events, zonal current anomalies were calculated for the years of interest. Figure 5-7a 
depicts the 13-years’ mean climatology of the zonal currents along the equator. Figures 
5.7b-d depict the anomaly of the zonal currents for the years with known positive Dipole 
Mode events. 






Figure 5-7. (a) Climatology of zonal currents (cm s'^) along the equator derived for 

satellite based observational data; (b) Anomaly of zonal currents during a 
positive Dipole Mode event year 1994, (c) Year 1997, and (d) Year 2006. 
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It can be clearly seen that all these years show consistent negative anomalies. This 
validates the previous studies with the observation that these jets slow down during a 
positive Dipole Mode event. This slowing event of the jets is analogous to the strong 
westward winds during a Positive Dipole Mode event. From the results, it was also seen 
that the intensity of negative anomaly is analogous to the strength of a particular Dipole 
Mode event. The biggest negative anomalies occurred during the latter part of 1997 
(Figure 5-7c) which is the strongest known positive Dipole Mode event in recent history. 
Negative anomalies can also be seen in the most recent positive Dipole Mode event 
(Figure 5-7d) which was in the latter part of 2006. The same can also be verified from the 
Dipole Mode Index (DMI). 
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Figure 5-8. (a) 13-year mean climatology of surface currents along the equator during 

November, (b) Surface currents anomaly during November 1997. 
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To see more closely how equatorial jets altogether reverse their direction of flow 
during a positive Dipole Mode event, vector plots for the equator region were generated 
showing the anomaly for the months with strongest DMI. Then, it was compared with the 
mean climatology of that particular month. For the Dipole event of 1997, the maximum 
DMI was recorded for November. Figure 5-8 illustrates the 13-year mean climatology of 
surface currents along the equator during November (Figure 5-8a) and, also, shows the 
anomaly during November 1997 (Figure 5-8b). It is dramatic to see how strongly the 
wind patterns during a positive Dipole Mode event affect the general surface circulation 
at the equator. 

F. EQUATORIAL JETS DURING NEGATIVE DIPOLE MODE EVENTS 

To analyze the behavior of the equatorial jets during known negative Dipole 
Mode events, zonal current anomalies were calculated for the years with known negative 
Dipole Mode events. Figure 5-9 illustrates the anomaly of the zonal currents for the years 
with known negative Dipole Mode events. From the figure, positive anomalies were 
evident in all three years. It reached a maximum in 1998 where the positive anomaly 
reaches up to 80 cm s'\ The results from satellite-based observational data confirmed the 
results of previous studies: it was observed that equatorial jets get stronger during a 
negative Dipole Mode event. The strong eastward jets may be caused by the strong 
eastward winds during a negative Dipole Mode event. Positive anomalies of the zonal 
currents can also be seen during the most recent year of negative Dipole Mode event 
during 2001 (Figure 5-9d). From the figure, it can also be observed that the intensity of 
positive anomaly is associated with the strength of a negative dipole event, which is 
represented by DMI. 

To further see the behavior of equatorial jets during a negative Dipole Mode 
event, vector plots were generated. During the negative Dipole Mode event of 1996, the 
strongest DMI was recorded during the month of October. Again, Figure 5-10 serves as a 
sample vector plot of what the equatorial jets look like during a negative Dipole Mode 
event. The figure contains the 13-year mean climatology of surface currents along the 
equator during the month of October (Figure 5-10a) and, also, illustrates the anomaly 
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during October 1996 (Figure 5-10b). From the figure, we can clearly observe positive 
anomaly during the month of October as still in the anomaly plot we can see that the 
currents are moving in eastward direction. 




Figure 5-9. (a) Climatology of zonal currents (cm s'^) along the equator derived for 

satellite based observational data; (b) Anomaly of zonal currents during 
Negative Dipole mode event year 1996, (c) year 1998, and (d) year 2001. 

57 











Mean Oct 



Longitude 

Surface Current Anomaly Oct - 1996 



45^ eo't 75*^ ao'fe 

Longitude 

Figure 5-10. (a) 13-year mean climatology of surface currents along the equator during 

October, (b) Surface currents anomaly during October 1996. 

From the above results and discussion, it is concluded that the satellite-based 
observational data can be used as an important tool to diagnose and even forecast the 
occurrence of an upcoming Dipole Mode event. The most useful thing about the OSCAR 
data set is that it is being updated on a regular basis (after every five days). Further, its 
availability is speedy compared to spot measurements. Although the data set may have 
certain inherent issues with the quality and accuracy, reasonably accurate results can still 
be obtained by applying certain mathematical techniques (such as OSD). Similar 
methodology, as we have utilized to analyze the Indian Ocean Dipole Mode events, can 
be used to investigate other large-scale ocean phenomena as well. 
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VI. PRACTICAL APPLICATIONS 


One of the most practical features of the OSCAR data is its near real-time and 
consistency (data available every five days). Because of this attribute of the data, its use 
for operational and tactical purpose becomes extremely useful: the results produced from 
this data set provide more realistic depiction of the synoptic surface currents in any area 
of interest. The data is equally useful for scientific purposes as well. In this subsection, 
we will briefly see how the knowledge of ocean surface currents can be practically used. 

A. OPERATIONAL USES 

In this subsection, there will be a brief discussion about the possible operational 
uses of studying ocean surface currents. 

1. Ocean Currents and Mine Warfare 

Mine warfare is an integral part of naval warfare. Naval mine is one of the most 
cost effective weapons as it is inexpensive, easy to conceal, easy to deploy, and requires 
almost no maintenance. These attributes of the naval mine make it the favorite for nations 
with small naval power. Mines can be used strategically and tactically for denying access 
and defense of vital facilities, such as ports, anchorages, and offshore structures. Mines 
can change the shape of a battle space and can effectively damage the operations of naval 
forces. One of the important types of mine is the drifting mine. Drifting mines are usually 
used to completely refute access of friend and foe to a battle space. The mines are easier 
to deploy but are very difficult to sweep and hunt. The synoptic ocean surface currents 
affect the surface mine drifting. Deploying locations of the mines could be decided by 
knowing the most recent surface circulation pattern in a particular region. Thus, better 
understanding of ocean surface currents will lead us to better and more effective mine 
warfare. 

2. Search and Rescue at Sea 

In both peacetime and battle space environments, search and rescue operations 
require meticulous planning and careful execution. Amongst other environmental factors, 
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ocean surface circulation is the most important factor in planning a search and rescue 
operation. Ocean circulation pattern could be used to estimate the likely position of men 
and material to be rescued. By doing this, rescue effort can concentrate in a smaller and 
more probable area of interception. The good knowledge of ocean surface circulations 
can tremendously reduce the search and rescue effort by channeling it to a smaller area of 
interest. 


3. Planning for an Economical Passage 

Surface ocean circulation patterns could be used to plan for a more economical 
passage. For longer voyages, without jeopardizing the shortest route, routes that are along 
the general flow of the ocean can be planned. By doing this, fuel consumption can be 
reduced and moves will be faster. 

4. Monitoring Chemical Pollutants and Oil Spills 

In an event of intentional or unintentional spills of chemical pollutants and oil, 
ocean surface currents could be used to determine the spill. This is especially true in the 
coastal regions with very fast surface currents, such as the Gulf Stream, Kroshio, and 
Somali Current. It tells us how fast to act in terms of remedial measures. Better 
understanding of ocean currents can also help in planning and executing such remedial 
measures. Thus, a good knowledge of ocean surface currents can not be undermined in 
context to monitoring unwanted liquid spillage in the ocean’s upper surface. 

B. SCIENTIFIC USES 

In this subsection, there will be a brief discussion about the possible scientific and 
technical uses of studying ocean surface currents. 

I. Ocean Currents as a Source of Energy 

In today’s industrialized world, the importance of sustainable energy sources 
cannot be undermined. The bigger economies of the world are always looking for newer 
ways to produce good amounts of energy. Ocean currents have the maximum strength 
near the surface and most of these are wind-driven. Power extraction from a flowing 
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source, such as ocean currents, depends on speed and density. Water is about 835 times 
denser than wind so, for equal area of intercepted flow, the energy contained in 12 mph 
would be the same as for a wind mass flowing at 110 mph (Technology white paper. May 
2006). Also explained in the Technology white paper is that the estimated world wide 
power in world’s ocean is about 5,000 GW with the power densities of about 15 KW/m . 
The Technology white paper. May 2006, also comments that if only 1/1000 of the energy 
in Gulf Stream is utilized, it will be about 21,000 times more energy than in the Niagara 
Falls. With this little comparison, it is clear that worldwide ocean currents house a 
tremendous amount of energy in them. Most of it is unused. 

2. Large-Scale Climate and Synoptic Diagnostics 

In terms of the OSCAR data set as the near real-time data of ocean surface 
currents as regularly available, with 1° x 1° resolution, it can, therefore, be used for large- 
scale synoptic diagnostics. Different large-scale features in the ocean surface currents 
could be analyzed regularly to see the variability. One of the examples in this context is 
the Indian Ocean Dipole phenomena. As seen in Chapter V, the occurrence of lOD can 
also be diagnostic via analysis of ocean surface currents. Besides this, variability in the 
large-scale eddies can also be monitored through ocean currents. Further, change in their 
general circulation behavior can be analyzed or used for large-scale climate diagnostics. 
In short, it can be said that overall ocean surface circulation constitutes an important 
factor in the overall climatology of a certain area. 
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VII. SUMMARY 


The OSCAR, from 1993 to 2006, satellite altimetry data, was used to analyze the 
seasonal behavior of the surface circulations in the Indian Ocean region. More focus was 
given to the three most dynamic regions of the Indian Ocean: Somali Current, North 
Arabian Sea, and Bay of Bengal. The same data set was also used to determine the link 
between the eastward flowing equatorial jets and the Indian Ocean Dipole Mode events. 
Zonal current anomalies along the equator during a Dipole Mode event showed that 
occurrence of such an event directly affects the flow pattern of these equatorial jets 
(Wyrtki jets). 

The raw altimetry (OSCAR) data was refined by using the Optimal Spectral 
Decomposition method. This study shows the tremendous ability of this method to 
improve the noisy OSCAR data with realistic depiction of the coastal boundary currents. 
The strength of this method can be judged from the fact that after reconstruction, the non- 
evident boundary currents, such as Gulf Stream, Kuroshio, and Somali current, were now 
evident and the flow patterns could be clearly seen. 

Using the reconstructed altimetry data, the mean seasonal variability of ocean 
surface circulations in the Indian Ocean region was studied. This study confirms, with the 
exception of a few small-scale features, the seasonal pattern of Indian Ocean surface 
circulations that was seen in previous studies. The reason why OSCAR data has a 
difficult time identifying small-scale features may be its resolution: it reduces down to 1° 
X 1°. The visuals produced from the data set show that, during summer monsoon, SEC 
and EACC supply the northward flowing Somali current. Part of the Somali current re¬ 
circulates a little south of the equator to form the “Southern Gyre;” the other part moves 
in the northward direction. Eurther north during this time, another gyre forms which is 
known as the “Great Whirl.” The Southwest Monsoon Current (SMC) south of Sri Eanka 
flows eastward. During this season, the East Indian Coastal Current (EICC), with a 
weaker strength, flows northward along the eastern coast of India, whereas the general 
flow inside the bay is eastward. Along the eastern coast of Madagascar, the South 
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Equatorial Current (SEC) splits up into two branches: one move northward as the 
Northeast Madagascar Current (NEMC); the other moves southward as the Southeast 
Madagascar Current (SEMC). 

During the winter monsoon, the East African Coast Current (EACC) meets the 
then southward flowing surface Somali current in a confluence zone. The two flows then 
supply the eastward flowing South Equatorial Countercurrent (SECC). At the eastern end 
of the SECC, the South Java Current flows southeastward. At this time of year, the 
southerly flowing Somali current gets input from the westward flow in the Arabian Sea. 
This westward flow is mostly contributed by the westward flowing Northeast Monsoon 
Current (NMC) south of Sri Eanka. After circulating around the Eaccadive High, the 
NMC also supplies the WICC. During the winter season, there seems to be a clockwise 
moving gyre in the Bay of Bengal. The EICC attached to it seems to have a southward 
flow. In Southern Hemisphere Region, there is a broad zonal eastward flow of the South 
Equatorial Current (SEC). This is driven by the Southeast Trades. The visuals produced 
from altimetry data show that ocean currents south of about 15° S do not show much 
seasonal variability; therefore, SEC shows almost the same flow pattern in both seasons. 

The behavior of the eastward-flowing equatorial jets (Wyrtki jets), in conjunction 
with the Indian Ocean Dipole events, was studied. The results showed that Dipole events 
directly affect the flow of the eastward equatorial jets. The unusual change in the flow of 
these jets is attributable to the changed flow pattern of the wind forcing as compared to 
their usual seasonal flow. This study confirms all previous studies: the eastward flowing 
equatorial jet weakens in a positive Dipole Mode event and strengthens in a negative 
Dipole Mode event. Both zonal current anomalies and vector anomaly plots depicted the 
same results. In short, it was concluded that satellite-based observational data is an 
important tool to diagnose and forecast the occurrence of an upcoming Dipole Mode 
event. 

The most practical feature of OSCAR data is its currency and consistent flow 

(data available after every five days). As a result, its use for operational and tactical 

purposes becomes extremely useful. This is because the results produced from this data 

set provide more current and realistic depiction of the overall surface currents’ 
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climatology in a particular area of interest. The data is equally useful for scientific 
purposes. For operational use, this ocean surface currents’ data can be used for mine 
warfare, search and rescue operations, planning economical passages, and monitoring 
chemical pollutants and oil spills. In scientific context, the data can be used for 
developing energy sources using ocean currents and for large scale climate diagnostics. 

For future research, this study recommends that OSCAR data can be used to 
analyze ocean surface circulations in other parts of the world as well. Although the 
satellite based data has inherent concerns pertaining to resolution and accuracy, it still can 
be useful in identifying mesoscale features in the ocean. With considerable refinement 
(OSD) and correct methodology this data can also be utilized to examine large scale 
climate phenomena, such as Indian Ocean Dipole events, El Nino, and La Nino. 
Comprehensive knowledge of the world ocean surface circulations will also serve vital to 
the Naval Forces in planning and executing of their missions over the sea. 
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